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A microscopic theory for equilibrium and non equilibrium 
relaxations in structural glasses is formulated. For all temper- 
atures below the glass transition the dynamics can be asymp- 
totically separated in a /3 - relaxation regime, which connects a 
quasi equilibrium with the aging regime and the aging regime. 
In the /3 - relaxation regime dynamic correlation are indepen- 
dent from spatial correlations. This is due to a self organized 
critical bifurcation scenario, which gives rise to one dynamic 
power law and a logarithmic law. It is shown that the fluc- 
tuation dissipation ratio in the aging regime is a universal 
property of the system independent from spatial scales. 

PACS numbers: 61.43. Fs, 61.20.Lc, 64.70. Pf 05.70.Ln 



I. INTRODUCTION 

Glasses are important solids, which are by definition in 
an out of equilibrium state. Nevertheless in many exper- 
iments they can be dealt with as if they were in a time 
independent quasi equilibrium state. E.g. measurements 
of static structure factors and dynamic neutron or light 
scattering functions give reproducible results seemingly 
independent of the waiting time 1^ after the production 
of the glass. On the other side it is known, that properties 
of glasses depend on the procedure with which they are 
obtained and that they change on very long time scales. 
The later phenomenon is known as aging of glasses and 
there exists a vast experimental literature on this non 
equilibrium behavior Q. So far it exists no microscopic 
theory of structural glasses, which contains the equilib- 
rium and the non equilibrium aspects of glasses, including 
the influence of their non ergodic character on the struc- 
ture and on the other side the influence of the structure 
on the dynamics, j^]. Although there are many theoreti- 
cal investigations on the aging behavior of simple models 
which are able to rationalize properties of simula- 
tions on models of simple liquids out of equilibrium |^-^ , 
they cannot, due to their simplicity, address structural 
aspects of the aging dynamics . For supercooled liquids 
in equilibrium a microscopic theory, the mode coupling 
theory (MCT) of the glass transition, exists, which suc- 
cessfully describes the dynamic of such systems at tem- 
peratures above the critical glass transition temperature 
1^,^. The applicability of this theory below the tran- 
sition temperature is not obvious, since it is based on 
equilibrium considerations, which are not fulfilled a priori 
for non ergodic systems ||l^. In this letter I will present 



a description of the dynamics of glasses, which unifies 
the non equilibrium aspects discussed in simple model 
systems and the detailed structural aspects of the (equi- 
librium) MCT of supercooled liquids. The theory allows 
to draw very general conclusions on the nature of aging 
in structural glasses. 

The goal is to develop a theory, which allows to cal- 
culate dynamic correlation and response functions in 
space and time for glasses. To describe the experimen- 
tal situation, only the initial equilibrium state is as- 
sumed to be known. The final glassy state is a re- 
sult of the dynamic reached for time t oo after the 
quench at t = 0. The variables relevant for the de- 
scription of structural relaxation are density fluctuations 
nq(t) = X^i 6xp(i(7r'j(i)) and longitudinal current fluctu- 
ations = X;j(#5)/(9w)exp(jgrj(t)). Here fj, pj 
and m are the position, the momentum and the mass of 
particle j. The density fluctuations contain the informa- 
tion about the spatio temporal distributions of the parti- 
cle, whereas the momentum fluctuations are necessary to 
capture possible effects of global momentum conservation 
and to describe the rate of the spatio temporal changes. 
The equations of motion are given by the underlying 
Hamiltonian dynamics and the coupling to a thermostat, 
which guarantees, that the chosen thermodynamic pa- 
rameters like density or temperature are kept constant 
after the quench. To be specific, I will choose in the fol- 
lowing a temperature quench from an initial equilibrium 
state with temperature Ti = 1/ {PikB) to a final tempera- 
ture T. The temperature can be measured with the help 
of the kinetic energy ksT = 2/{3N){J2jP'j/i'2m)), where 
N is the number of particles. The average (. . .) is an 
ensemble average over the initial non equilibrium ensem- 
ble pne = exp(-/3E/p')/(2m) - A E.<, Viln - r,\)), 
where V{\ri — fj\) is the interaction potential. This 
ensemble guarantees, that the value for the temper- 
ature, measured with the help of the velocity distri- 
bution is Ub/P and that the spatial structure is the 
same as that of the liquid at the initial temperature 
ks/Pi- The dynamics is given by a time dependent 
differential operator L{t). For a variable A{fi,pi) the 
equation of motion is: dA/dt — iU{t,Q)L{t)A{0) := 

U{t, 0) {iLoA - a(t)p,{0)dA{n,pi)/dp,{0)) where 

C/(t, 0) — T+exp(i j^L{t')dt') is a time ordered prod- 
uct (time increases from left to right). Lq = —i{H, .} is 
the Liouville operator, and {., .} is the Poisson bracket. 
The Gaussian thermostat a{t) is adjusted such that 
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(Xljbj (i))/(3iVm)) = 1//3 for all times t. Given these 
equations of motion for the phase space variables, one 
has to derive the equations of motion for the experi- 
mentally relevant correlation functions and, since a pri- 
ori the validity of the FDT is not guaranteed, also for 
the susceptibilities. A suitable method to achieve this, 
is a non equilibrium projection operator formalism, us- 
ing nq{t) and jq{t) as slow variables |13|. To obtain 



an exact equation for the susceptibilities, it is necessary 
to formulate the projection operator in terms of Pois- 
son brackets instead of correlation functions With 
such a projection operator and the specified iso-kinetic 
dynamics the equation for susceptibility Xq {t-w , t) = 
eit - t^)[ng{t^)\ng{t)]/N := {{n*{t^),ng{t)})/N and 
the correlation function Cq{tw,t) = {nq{tw)\nq{t))/N :— 
{'^*q{tw)\nq{t)) / N ^ where t^^ is the waiting time after the 
quench, can be derived as 

^Xq{tw.t) = ^l-Xq{t^,t)[iML{i)3q{t)] (1) 

-ima{t)—Xq{tu],t)+m J dt'xqitio,t')Y.q{t' ,t). 

The memory function Y,q{t' ,t) is the average of the 
Poisson bracket of the true force dX t = t' and the ran- 
dom force Fq{t',t) between t', and t. This random force 
describes the dissipation of the gradient of the longitu- 
dinal stress at t' in non-hydrodynamic modes between 
t' and t. The correlation function Cq(tw,t) depends on 
all times between and t, since the time translational 
invariance is broken. The equation for it reads 

0'^ Tfl 

— Cq{t^,t) = -l-Cq[t^,t)[h{i)\L{t)3q{i)\ 

d /■* 

- nna{t) — Cq{t^,t) + m j dt' Cq{t^,t')Y.q{t' ,t) 
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The last two lines in Eq. ^ describe the correlation 
to the initial state. The memory function Mq{t^,t) is 
an autocorrelation function of random forces Fq{Q,t.^) 
and Fq{Q,t). The equation for the time dependent equal 
time correlation function Sq{t) = {nq{t)\nq{t))/N , which 
describes the change in the structure of the liquid after 
the quench, is easily derived. 
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2^5,(i) = {jq{t)\jq{t))/N - —Cqit^,t)\t^^f (3) 

If we make the assumption, that the velocities of 
different particles are as in equilibrium uncorrelated 



at every time t, equation simplifies further, since 
ijq(.'t)\jq(.'t))/N = '^/(Pm) in this case. 

The set of equation are formally exact. If the 

initial state is an equilibrium state i.e. if there is no tem- 
perature quench, it can be exactly shown [Q, that the 
correlation function and susceptibilities fulfill the FDT 

Xqitw,t) = Xq'it - t^) = -{l/T)dCl'i{t - t^)/dt (4) 

and the memory functions I]q(i„,,t) and Mq{tw,t) obey 
a FDT of the second kind 

i) = - t^) = -{l/T)dM'q'i{t - t^)/dt (5) 

The three equations (|^, ^) and are reducing then 
to one equation for the correlation function ||T^, which 
was the basis for the equilibrium MCT [ p^ . 

A solvable theory is obtained by approximating the two 
unknown functions T,q{tw, t) and Mq{tw,t). The simplest 
non trivial approximation consistent with the mathemat- 
ical requirement for the memory functions is a one 
loop type or mode coupling approximation. It has the 
form 



k,p 

Cqkp{{Sq{tu])]) Cqkp{{Sq{t)])Xk{tw,t)Cp{tw,t) 

Mq{t^,t) = <zV(2™2^2) J2 6. .^z 

k,p 

Cqkp asqit^..)}) Cqkp 

{{Sq{t)})Ck{t^,t)Cp{t^,t) 



(6) 



(7) 



The vertices Cqkp are functionals of the one time corre- 
lation function Sq(t). They are symmetric with respect 
to exchange of p and k. Under this general condition, the 
FDT of the second kind (^) is automatically fulfilled in 
equilibrium. The explicit form of Cqkp{t) will determine 
the theory for the structure factor Sq{t). For asymp- 
totically long times i^, , i — > cxd all correlation functions, 
which depend only on a single time can be replaced by 
their asymptotic value. If the quench was performed from 
an equilibrium liquid state, the functions Cq(t^,0) and 
(rig(itu)ljq(O)) in Eq. (||) will vanish for tw ^ oo. It can 
also be shown [Q, that the influence of the thermostat 
a{t) vanishes at least as 1/t and therefore the equations 
for the correlation functions and the susceptibilities in 
the long time limit depend only on the asymptotic value 
of the structure factor S'^ and the vertices Cqkp{oo). 

If there were no wave vector dependencies, the struc- 
ture of the resulting equations were similar to the equa- 
tions for the spherical p - spin model [|16|,| in the same 
asymptotic limit. As in this model I will make the ansatz, 
that there exist below the glass transition two differ- 
ent dynamic regimes. The FDT regime is characterized 



by the hmit Cf*^^ (r) 



lim* 



> Cq{tw,t^ + r) with 



(-l/r)dC^^^(r)/dT. The aging regime is 
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defined by the limit Cq{X) = limt^oo Cq{Xt, t). This scal- 
ing behavior is called simple aging. There are more gen- 
eral scenarios which could so far only be ruled out for 
the p - spin model For notational simplicity I will 
concentrate on the ansatz of simple aging. I also will 
make the ansatz of a generalized FDT ansatz, known to 
be correct for the p - spin model It assumes, that 
the susceptibility Xg(A) = limt_»oo tXq(Ai, is related to 
the correlation function by Xg(-^) = {Xq/T)dCq{X)/dX 
psf . Astonishingly it turns out that only a wave 
vector independent fluctuation dissipation ratio (FDR) 
< X{T) < 1 will give compatible equations for correla- 
tion functions and susceptibilities in the aging regime. 

With the definitions of normalized quanti- 
ties ^q{t^,t) = Cq{t^,t)/S^, Nq{t^,t) = Xq{t^,t)/S^, 

fq - Fq/S^, where Fq = lim.^oo ^^(r) - (7,(1) and 
Vqkp = '5^,fc+j?c^fcp(oo)S'^S'jf', the following equations for 
the asymptotic structure factor is obtained 

1/5,°° - m/34 - l/2Y^Vqkp { 1 - fqfkfpil - X)} (8) 

k,p 



/./(I - A) = E "«fep/fc/f (1 - (1 - - A)) (11) 

Since it can be shown, that X — 1 for T = only, Eq. 
( |rT| ) reduces to the equation of the MCT of supercooled 
liquids only at T = Tc. For T <Tc\t is modified due to 
the breaking of the FDT {X^l). 

A deeper insight in the meaning of the wave vector 
independence of X can be obtained by investigating the 
behavior of(j)q{tw, t) for times, where it is close to the non 
ergodicity parameter Fg. A detailed analysis shows , 
that the Eqs. (^) and (^0|) can be simultaneously solved 
by the ansatz 4>q{tw,t) = fq + eqG(tw,t) up to quadratic 
order in G{tw,t). Here is an eigenvector of a matrix 

Dqk 



= J2Dqk'G'k,il) with 



(12) 



Here 



Unit ^oo 'irn[jq{t)\LQjq{t)]/q'^ is a bare wave 
vector dependent velocity, which will be renormalized 
by contributions from the vertices Vqkp- The equa- 
tion for the correlation function in the FDT regime i.e. 



-/3 d<j,^D^{T)/dT = xr'^(^) is given by 
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2m/3 4^ 

kp 



qkp 



Dqk' = - (gV2™/3) E"9fep(l - fkfp))Kk' 

kp 

+ iqy2mP)vqk'pfpil - fq). (13) 

Eq. ( [T2I ) is the solvability condition for Nq{X) — 
{X/T)dGq{X)/dX by setting A = 1 in Eq. (0). 

The equations for G{tw,t) in the FDT and the aging 
regime follow by multiplying the resulting equations from 
left with the left eigenvector e, to eigenvalue zero of Dqk' 
i-e- J2„^qDqk' = 0: 



(9) 



= A (G^OI^)2(t) - f dT'G^''^{T-T')G^^^{T') 

dT Jo 



2m/3 



k.p 




,rf0r^(r'), 



dr' 



ir^x<^^)(r-r') 



+ fqfkfpil ~X)} 

With Eq. (||) and (||) it follows immediately that Eq. 
(^ reduces to the standard MCT equation for the su- 
percooled liquid above the glass transition temperature 
Tc, where fq = and Xq — 0. It is interesting to note, 
that like in the approach of Zaccarelli et al. |l^, the 
mode coupling approach contains implicitly a theory for 
the structure factor of liquids. Below the glass transition 
the structure of the glass depends on the non ergodicity 
parameter i.e. the non ergodic nature of the glass also 
influences the equal time correlation function. 



0= (A/X)G2(A) 



ds 



Gis) 



where 



A = 



(1 - /<z)efeep 



'^qkp mfj'^^qkpe-qekfp 



(14) 



(15) 



(16) 



= + 



k,p 



The microscopic structure of the system only enters via 
the exponent parameter A, which is a complicated func- 
tional of the asymptotic structure factor S^. This be- 
havior is very reminiscent of the (3 - relaxation phenom- 
, ^ ena in supercooled liquids above the glass transition. Eq. 

l/(2m/3) Vqkp |(1 - fkfp))Nq{X)) (10) (|l|) has exactly the form of the P - scaling equation of 

the MCT of supercooled liquids at the critical tempera- 
ture ||2^]. In this theory the glass transition is a critical 
bifurcation of codimension 1 , in which the dynamics close 
to the transition temperature and for times, at which 
the correlation function is close to the non ergodicity 
parameter, is completely determined by the critical di- 
rection of the bifurcation. In the theory presented here, 
the additional solvability condition for aging Eq. ( p^ ) 



+ iVfc(A)<^,(A)(l - fq) + 1/13 r ^Nq((-))Nk{s)4>p{s) 

Jx s s 

The Eq. for the non ergodicity parameter can be ob- 
tained from the t 00 limit of (||) and it is identical 
with the initial value 4'q{l)- With the help of Eq. (^) it 
can be written 
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forces the system always to be at the critical bifurca- 
tion for temperatures below T^. Since this bifurcation 
is still of codimension 1, there is only one temperature 
dependent number to guarantee this self organized crit- 
icality, namely the fluctuation dissipation ratio X{T). 
This is the deeper reason, why X does not depend on 
wave vectors. Within the presented theory, aging is due 
to a self organized critical bifurcation scenario. Because 
of this criticality Eq. (Q) and (15) can be solved ex- 
actly for all temperatures below the transition tempera- 
ture G^^-^(t) = T~°- for t* T To where tq is some 
microscopic correlation scale and t*(ttu) is the time where 
4>q{tw, t) = fq [^3| . The exponent a is determined by A: 



A(r) = r2(l-a)/r(l-2a). 



(17) 



where F is the Euler Gamma function. The form of this 
equation is known from the /3 - scaling theory of super- 
cooled liquids (2^ . But contrary to this theory Eq. ( |l7| ) 
holds for all temperatures below the glass transition with 
a temperature dependent exponent parameter A and thus 
a temperature dependent exponent a. Equation (^5|) is 
solved by G{\) = -(-logA)'' for A < 1. Here ta{tyj) is 
a time characterizing the final decay of the correlation 
function The exponent b is determined by an equation 
similar to ( p7| ) 



A(T) = X{T) T^il + b)/T{l + 2b) 



(18) 



At the transition Eq. ([ij) reduces to the Equation for 
the von Schwcidlcr exponent of the (3 scaling theory . 
Below the transition relation (^8|) was also found in mean 
field theories |^l|, in which the dynamics reduces for all 
times to a single component theory. For models, in which 
simple aging does not hold, A — t^/t has to be replaced 
by A' = h{tw)/h{t), where h{t) is a yet unspecified func- 
tion (see e.g. the discussion in [0). It is important to 
notice, that the decoupling of wave vector dependence 
and time dependence in the present theory does not ap- 
ply once the condition &{X} ^ G'^(A) is violated. In this 
case wave vector dependencies of the correlation func- 
tions are to be expected and the solution of (|l^) can only 
be found numerically {X(T) of course remains wave vec- 
tor independent). The predicted factorization property 
4>q{tw,t) = fq + eqG{tw,t) was indeed found in simula- 



tions ||7|. The presented theory allows for a systematic 
theoretical investigation of the dynamics of the glassy 
state of simple liquids as e.g. the colloidal hard sphere 
system, in which indication for aging was already ob- 
served ||2^] or the Lennard Jones systems of 1^,^. Es- 
pecially it the foundation for a theory of the structure 
factor of glasses by specifying the vertices Cqkp (t) in Eqs. 
(^ and (|^). This will be part of future work. 
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